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Abstract 

To check the consistency of positivity requirements for the two-point 
correlation function of the topological charge density, which were iden- 
tified in a previous paper, we are computing perturbatively this two- 
point correlation function in the two-dimensional 0(3) model. We 
find that at the one-loop level these requirements are fulfilled. 



1 Introduction 

In a preceding paper pQ we analyzed the implications of physical positivity 
('unitarity') and positivity of the topological susceptibility for the two-point 
correlation function of the topological density. In particular we found that the 
short distance singularity has to be softer than tree level perturbation theory 
would predict. In this article we verify for the two-dimensional 0(3) nonlinear 
sigma model that one-loop perturbation theory is indeed in agreement with 
this requirement. Such a softening of the short distance singularity is already 
predicted by including the running of the coupling constant in tree level 
perturbation theory E] . But since we are dealing here with a composite 
operator, this check of internal consistency of perturbation theory is quite 
nontrivial. 



1 



2 Topological charge density correlator 



The action of the (ferromagnetic) 0(3) spin model on a 2-dimensional square 
lattice A C Z 2 is 

5=-^^-Sy, (1) 

where the sum is over nearest neighbor pairs (XY) of lattice sites. The 
partition function reads 




with measure 

[DS] = Hd*S x 5(S*-l), (3) 

X 

imposing the unit norm constraint for all spins, i.e. S G S 2 C M 3 VX G A. 

We will consider the thermodynamic limit as the L —>■ oo limit of the 
model on an L x L square lattice with periodic boundary conditions, so that 
translation invariance holds. The origin is then located at the center of A. 
Unit lattice vectors are called 1 and 2. 



2.1 Definitions of the topological charge 

We want to study the two-point function of the topological charge density 
qx, to be denoted by £ x (with x =0X), 

£x = Go,x = (qxqo) = \ J [BS}q x q e- s . (4) 

The topological susceptibility is defined by 

Xt = ^Go,x- (5) 

x 

In this paper we will consider two lattice definitions of qx'- a 'field theo- 
retical' one and the geometric one due to Berg and Luscher. Both make use 
of a triangulation of the square lattice Z 2 as in figure ^ I n both cases, the 
topological charge density is written as 

qx = q x + q\, (6) 
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where 

& = f(X,X + l,X + 2), 

(7) 

q b x = f(X + l + 2,X + 2,X + l), 

with / : A 3 — > K a certain function. 

In P] we used a symmetrization of this definition by introducing in addi- 
tion a 'mirror' triangulation 

q x = f(X,X + l,X + l + 2), 

(8) 

q d x = f(X,X + l + 2,X + 2) 

and defining q x = oi^x + 9x + <?x + Qx)'i this was necessary there in or- 
der to maintain reflection symmetry. Here, however, this symmetrization is 
unnecessary, because it has no effect in perturbation theory. 

The 2-point function Q x = £ ,x thus has the form 

& = (qxqo) = <(<& + QxM + ?o)> = E S * ^ 

i,j=a,b 

with 

G2 = {<tx<&), i,j = a,b. (10) 
To make the notation more transparent, call 
(Z u Z 2 , Z 3 ) = (0, + 1, + 2), 

(W 1 ,W 2 ,W 3 ) = (X,X + l,X + 2), 

(Z[, Z' 2 , Z' 3 ) = (0 + 1 + 2,0 + 2,0 + 1), 

(W{, W' 2 , Wi) = (X + 1 + 2, X + 2, X + 1) 



11) 



and denote Z, = SV^ and analogously for W, Z', W (see figure [TJ). The 
indices just introduced are denoted by i, j, k, etc. 

In the following sections, only formulae for q\ and Q® a will be presented. 
The full topological charge density can be recovered by means of (0) and (J7J), 
and the full two-point function by adding the other three contributions in © , 
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X W x W 2 X+l 



0+2 Z' 




Zi Z 2 0+1 
Figure 1: Triangulation of the lattice plane. 



obtained by substituting Z'W, ZW and Z'W for ZW, respectively. More- 
over, a numerical factor will be introduced in the definition of the correlator 
so as to work with simpler expressions: 

£c = (8tt) 2 £ x . (12) 



2.1.1 Field theoretical definition 



We first consider a symmetrized version of the field theoretical (FT) definition 
jl] of the topological charge density 



FT, a 



— Sx ■ (Sx+i A Sx+2), 



which gives rise to the two-point function 



'FT,aa 



Z x ■ {Z 2 A Z 3 ) W x ■ (W 2 A W- 



( det(^ • Wj] 
\ hi 



(13) 



(14) 



The perturbative treatment of the problem begins with the 0(3) — > 0(2). 
decomposition of the spins, 



Substituting in the determinant, we get 
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^ FT - = ^^et(.1.< 



77"^ ^ijk^tmn (zi- Zj • W m Z k • W r 

ijktmn 



+ 77o9 ^ijk^imn (a/1 7T A / 1 % ' "'m ^ " W„ 

2 ^ yS. \ V P V ^ 

(16) 

The perturbative expansion of ()16)1 consists of a Taylor expansion in 
Only integral powers of contribute. For our purposes, it will be enough 
to keep terms up to and including order /3~ 3 . Then 



(0) 

23 

ijktmn 

(1) 



ijkimn 



Zj ■ w m z k ■ w n 



1 _ _ _ _ _ _ \(°) 
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- - ( (3zi 2 - 2zi ■ w e + 3w e 2 ) Zj ■ w m z k ■ w 



+o(r 4 ), 

(17) 

where, in general, perturbative contributions to magnitude M. are denoted 
by 

M pcrt - — M {n) . (18) 
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2.1.2 Berg and Liischer's definition 



Berg and Liischer's (BL) definition of the topological charge density was 
introduced in jH]. With each elementary triangle (X, Y, Z) of the chosen 
triangulation one associates a multiple of the signed area of the minimal 
spherical triangle determined by (Sx, Sy, Sz) on the unit sphere. Explicitly, 

BL,a 1-1 S X ■ (S X +1 A Sx+2) / 1n x 
q x — 7T- tan =; - - =; =; — . (19) 

^ 1 + Sx ■ Sx+i + Sx+i • Sx+2 + Sx+2 • Sx 

The corresponding perturbative two-point function is, up to and including 
order /3~ 3 , 

PvBL,aa;pert. 

yx 

[1 + Zx ■ z 2 + z 2 ■ z 3 + z 3 • + m ■ W 2 + W 2 ■ W 3 + W 3 ■ W x ) 
+ o(/?- 4 )- 

(20) 

Substituting (JT5|) in (j2D|) and Taylor expanding the result, we arrive at 

£^BL,aa;pert. ^FT,aa;pert. _|_ ^Q aa (21) 

where 

1 / / - ■ ■ (0) 

ijkimn \ \u=l u<v 

+ o(/r 4 ). 

(22) 

will be referred to as 'BL difference.' 

2.2 Perturbative computation of Q x 

In this section we use the BL topological charge density (fTUj) and will com- 
ment on the difference with the FT definition when appropriate. Therefore 
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we drop the BL labels and write 

Gr l - = j 2 Gl + j^ + 0{(3-% (23) 

For bookkeeping purposes we do the perturbative computation of the 
density correlator out keeping factors N — l explicit, as for an O(N) problem, 
although the topological density has a meaning only for N = 3. 

The only vertices relevant for Feynman diagrams up to and including 
order (3~ 3 are order f3~ l : 

• Two-7r vertex: 

V& = i (l - ^y 1 ) 5 k % w . (24) 

• Four- 7? vertex: 

-n z S Z wTu{s ke S mn + 5 km 5 en + 5 kn 5^ 

+ 1 d ZT ozwotuo o +o zw d Z TOwuO o +o zw o Z uOwtO I , 

(25) 



\rklrnn 

V ZWTU — o 



where subindices are lattice points, superindices are 0(N) indices, nz = 4 is 
the number of nearest neighbors of point Z, two-index deltas are Kronecker, 

[I if Z = W = T = U, ( . 

dzwru - otherwise (26) 



and 



,, i 1 if Z, W nearest neighbors, , . 

° zw ~ i otherwise. { ' 



The explicitly volume-dependent factor in (J24)) comes from fixing one 
spin using the Faddeev-Popov method as in [0] to eliminate the divergent 
zero mode that would otherwise appear in the Gaussian measure defining 
tree level perturbation theory with periodic boundary conditions. 
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2.2.1 Green's functions 



In order to compute the Green's functions relevant to the correlators being 
studied, we start from the corresponding quantity on the Lx L square lattice 

and put x =0X= (x, y), —L/2 < x,y < L/2. 

r 0,2 V' cos(qp-x) 

°' X 2L 2 L> 2-cos(ap x ) -cos(ap y ) ' 1 ' 

where the momenta 

2tt 

P = (Px, p y ) = ~j- {n x , n y ), n x ,n y eZ (29) 

are summed over the first Brillouin zone — L/2 < n x ,n y < L/2, the prime 
meaning exclusion of the zero mode p = 0. 

The thermodynamic (L — > oo) limit of (|28p does not exist, since the sum 
diverges logarithmically in L. The divergence is independent of x, however, 
and the thermodynamic limit of the difference G 0t x — C 0i x — G does exist. 
According to a general argument due to David [7] , the perturbation expansion 
of invariant observables should be free of infrared divergences; the first sign 
of this is the fact that everything is expressed in terms of Go,x, which has a 
well defined thermodynamic limit given by the integral 

r -r r 1 f + \ [ + \ cos(q ■ x) - 1 

G , x = G , x -G = — — - / dq x / dq y ; 30 

2(27r) J J_ n J_ n 2 - cos q x - cos q y 

(another possible source of infrared divergence is the sum over intermediate 
lattice sites, which will be addressed in section 12 ,4j) . Obviously G = 0, 

This subtracted propagator was studied by other methods in [Hj, which 
we will use as a check for our results. 

Due to parity, the numerator in the integrand of (jHOjl can be replaced by 
cos(xq x ) cos(yq y ) — 1. Integration in q x can be performed by the calculus of 
residues, with the result 



1 f +1 dt (2-t + V3-& + t 2 y x T y (t) - 1 



G °' X ~ 27 U VT=? V 3-4t + f • (31) 
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Here T y (t) is Chebyshev's polynomial of the y-th kind. 

It is straightforward to check that (|3Tj) behaves asymptotically as 

Go,x = - ^ In |x| - c - A + 0(|x|- 2 ) (32) 



d = 0, (33) 



with c, d constants. According to jH], 

_ + 3 In 2 

where is Euler's constant. 
Similarly it can be seen that 

Go,x+v + GV-v-2G ,x = - f-A- + Odxl-') (31) 

7T \2X^ X 



and 

G G +0 G - 1 v 2 -w 2 -2(vx) 2 + 2(w.xf 

2ix |x| 4 

(35) 

for large |x| and |v|, |w| = 0(1). To see this, one subtracts from the Fourier 
representation of the lhs of (J34j) and ESI the corresponding continuum expres- 
sion with a suitable smooth high momentum cutoff; this difference, being 
the Fourier transform of a smooth (C°°) function will decay faster than any 
power of |x|. The cutoff on the continuum expression, on the other hand, 
does not affect the leading long distance behavior, which is obtained straight- 
forwardly. Expressions (J34|) and (J35|) . in which the constants c and d in (|32j) 
do not appear, will be useful in the computation of the charge correlator. 

In the following, we will denote Green's functions by Go,x even when 
subtracted. 



2.2.2 Tree level 

The tree level contribution to is given by 

1 v-^ / \ (0) 

^x' aa = 9 £ijk£emn(Zj ' w m z k ■ w n j , (36) 



ijkimn 
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Figure 2: Tree level diagram structure. 



the correlator in the rhs being evaluated in the Gaussian measure with one 
spin fixed. The corresponding diagrams have the structure of figure (J2J). 

Contracting spin indices, 



ijklmn 

+ 2 ~~ 1) £ ijk £ lmnG ZjZfiw, 

ijklmn 



(37) 



ijklmn 



The 3x3 matrices Gz.z. an d G\y.w. are symmetric, and therefore the 
second term in the rhs vanishes. Playing with indices, we can recast equation 
(13 7|) in the form 

gT,aa = ![(#- I) 2 - (JV - 1)] J2 e l3 kSlmnG Zj W m Gz k W n , (38) 

ijklmn 

This can be written in terms of differences G v — G w with |v — w| < y/2 
(which restricts IR divergences to summation effects, the IR divergences of 
Green's functions being cancelled in the differences): 

Ql aa = ~\[(N- I) 2 ~(N- 1)] 

(39) 

x 2^ £ijk£imn(G Zj w m - G Zk w m )(G Zk w m - G Z k W n ) 

ijklmn 

We can now perform the sums over all indices: in general, for a 3 x 3 
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matrix a^, 

^2 £ ijk^imnajma kn = ^ { [( tr a ) 2 - tr ( fl2 )] S M + 2(a 2 )^ - 2(tr a)a&} 

ijklmn it 

= 3 [(tr a) 2 - tr (a 2 )] + 2 s(a 2 ) - 2(tr a) s(a), 

(40) 

where s(a) denotes the sum of all entries in a. 

This identity can be obtained as follows: one starts from the well-known 
identity 

mdet(l + zA) = tr m(l + zA) (41) 
and expands it to order z 3 to obtain 

trA A A A A = -tiA 3 - -tr A tr A 2 + -(tr A) 3 . (42) 
3 2 6 V ' y ' 

By 'polarization', i.e. the replacement of A by xA + + zC on both sides 
and comparing the coefficient of xyz Eq. (}4T)|) follows if we put B = C = a 
and take for A the matrix with all entries equal to 1. 

In our case, 

Ql aa = -\ [(N - l) 2 - (N - 1)] 

x {3 [(tr G z . w f - tr (G| w )] + 2 s(G| w ) - 2(tr G z .w. ) s(GW. ) } , 

(43) 

2.2.3 One loop 

The one-loop (order /3~ 3 ) correction to the correlator can be decomposed as 

gL,aa = ^ya + ^ya + ^jj^aa + (jyya + (yjaa^ (44) 

according to the structures depicted in figure ©• 

• (I) contains the contributions of the 4-legged, disconnected graphs with 
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Figure 3: One-loop diagram structures. 



one 2-vertex: 

ar = [(N-if-(N-i)} 

(45) 

X £ i3k £ imnGz j W m ^ G Zk pG Wn P- 

ijklmn P 

• (II) contains the X-shaped contributions: 

(Iir = -\[{N- if — (N — 1)] £ ^rnn 

ijklmn 

x ^2 { G z 3 pGw m p - G Zj QG WmQ ) (G Zk pG WnP - G Zk QG WnQ ) . 
(P,Q) 

(46) 

• (III) gathers the contributions of 4-legged graphs with tadpoles: 
(HI)- = _ EJZI [( N - if - (TV - 1)] £ ijkStmnG Zj W m 

ijkimn 

x ( Gpp ~ G QQ) ( G z k pG WnP - G ZkQ G WnQ ) 
(P,Q) 

- [(N - I f — (N — 1)] £ijk£emnG ZjWrn ^ 

ijkimn 

x z^Z \GppG Zk pG Wn p + G QQ G Zk QG Wn Q 
(P,Q) 

-Gpq (G ZkP G Wn Q + G ZkQ G Wn p) ]■ 
The first term of (III) vanishes by translation invariance. 
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• (IV) consists of the six-legged graphs contributing to the correlator 
according to the FT definition of the charge density: 

N — 3 

(IV)- = — _ [(N - l) 2 — (N — 1)] detGW. 
N-l 



[(N - l) 2 - (N - 1)] ^ £ ijk£emnGooG Zj W m Gz k W n 

ijk£mn 



- [(iV-l) 2 -(iV-l)] 



X ^2 £ ijk £ imnG Zj W m (G Z,Z k G ' Zi W n + G Zk W e Gw e W n ) 
ijklmn 

The determinant term in (IV) aa has a vanishing prefactor. 

• (V) is the BL difference, also consisting of six-legged graphs: 



(48) 



(V) aa = - ^ [(N - l) 2 ~(N- 1)] 



y~] (G Zu z v - G 00 ) 



x 



^ £ ijk £ e mn Gz j w m Gz k w n 



ijklmn 



I [(TV - l) 2 - (TV - 1)] ZiJkeemnGz, 



ijklmn 



x l Gz ^ z k (Gz v w n - G Zu w n ) + G Zk w u (Gw n w v - G Wn w u ) 



u,v 



(49) 

As happened with the tree level expression, the one-loop correction can 
be rewritten in terms of differences G v — G w with |v — w| < y/2. This works 
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independently for (I), (II) and (V): 



(i-?LS\[(N-ir-(N-D] y, « 

ijklmn 



ijk^lmn 



■{G Zj w m - G Zk w m ) ^2(G Zk p - G Zl p){G Wn p - G Wm p), 



(50) 



(n,r = ^ [(JV - 1) 2 - (N - i)] * 



ijk^l inn 



ijkimn 



x \^ z i p ~ GzjQ) - {GwraP - G Wm Q)][{G z k p - G Zk o) - (G WnP - G WnQ )} 
(p,Q) 

x [(GzjP - G Zk p) + (G Zj Q - G Zk o) - (G Wm p - G Wn p) - {G WmQ - G WnQ )f, 

(51) 



(V d )" = [(N - I) 2 — (TV — 1)] 



^ (Gz u z v - Gqo) 



x ^2 £ ijk£emn(G Z] w m - Gz k w m ){Gz k w m - Gz k w n ) 



ijkimn 



I [(TV - l)^ _ (TV - l)] e tJ keemn(G Z] w m -Gz kWm ) 



(52) 



ijkimn 



X 



+ {G Zk w u - G Zk w v )(G WnWv - G WnWu )\, 

while (III) aa and (IV) aa cannot be brought to that form. However, the sum 
(III) aa + (IV) aa can, because, using translation invariance and the fact that 
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(III d ) aa - 2G o ^2 £ ijk £ emn(G Z:j w m - G Zk w m ){Gz k w m - G ZjWn ) 



ijklmn 



2(Gqi — Goo) £ijk£imn{GzjW m — Gz k w m ) 



ijktmn 



7, {Gz k p - Gz k Q){G Wn p - G w „q) 



(p,Q) 



(53) 



g ^ 

+ g(^01 ~~ ^Oo) /J £ijk£emn(G Zj W m ~ G Z k W m ) 
ijklmn 

x (Gz k p - G z ,p){G WnP - G WmP ) 
(p,Q) 

and 

(IV d ) aa = -2G o 2J ^'^(GzjWm - Gz fc w m )(G ? z fc Ty m - G Zj w n ) 

ijklmn 

~ ^2 £ ijk £ lmn(G Zl Z k ~ Goo)(G Zj W m ~ G Zj W n ) (G Z,W n ~ G Zk w„) 
ijklmn 

~ ^2 £ ijk £ lmn(G WeWn - G 00 )(G Zj W m ~ G Zk W m ) (G Zj W e ~ G Z k W n 
ijklmn 

(54) 

The first terms in the rhs of (|53|) and (|54j) cancel. 
2.3 ( a b notation 

Green's functions appearing (up to 1-loop order) in the perturbative expan- 
sion of Q x are of the form G v +w, where 

• v is one of ( =0X, £ =P0 (where P is summed over all lattice points), 
rj =PX and k = =00, 
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• w is a linear combination nil + 7^2 of unit lattice vectors with n^ 
integers (\rii\ < 2 up to 1-loop order.) 

The following notation proves useful when applying computer techniques 
to this perturbative problem: 

G*v+w > ^word- (55) 

Here v is one of the letters (, £, r] or k according to the vector v. The 
subscript is the minimal word of the form a na c nc b nb d nd , rti > 0, matching 
vector w under the rule 

a na c n cb n bd n d _^ ^ _ n< ,) 1 + _ (56) 

(obviously, only words of the forms a na b nb , a na d nd , c nc b nb , c nc d ni can be min- 
imal.) For instance, 

Co,X+l+2 ~~ > Cab, Gp^ 0+1-2-2 — ¥ £,add, G^O-l — ¥ K c- (57) 

This notation is way more readable than the standard notation, and al- 
lows us to cope with the full form of the density correlator, including all 
four terms in Q. AC program was used to generate the full analytic forms 
in Cab notation for each contribution (tree level and one-loop contributions 
I through V), its output being simplified by means of a computer algebra 
program. The simplification thus achieved is dramatic, witness the results in 
next section. 



2.3.1 The correlator in ( ab notation 



The full tree level 2-point topological charge correlator between and X, 
both in standard and Cab notation, reads: 



— 2 (Go,x+i+2 + G*o,x-i-2 — 2Go,x)(G'o,x+i-2 + Co,x-i+2 — 2Gq, 



x. 



— (Gti,x+i — Co,X+2 + Co,x-i — Ccx-2) 



(Cab + Ccd ~ X){Cad + Ccb ' 2C) " (C« -Cb + Cc-Cd 



(58) 
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From this expression and the asymptotics (J34|) and (J35|) of the differences 
of Green's functions, the leading term of the asymptotic behavior for can 
be easily computed: 

S = -^ + 0(|xr 5 ). (59) 

The full 1-loop contributions will only be presented in Cab notation (can- 
celling terms in (III) and (IV) are omitted): 

• the graphs with 2-vertices: 
(W = 2(l-|) 

(Ca - Cb + Cc ~ Cd) - &)(Vab ~V) + (U - 0(Va - Vb)] 

p 

- (Cab + Ccd ~ 2C) ^(£a - &)(*?„ " Vb) - (Cad + Ccb ~ 2{) J^(U - 0(Vab 
P P 

(60) 

• the X-shaped graphs: 

(IIC) = 2 - 0(U - U) - (£« - 6)& ~ Caab)} 

P 

X [(Vab ~ V)(Vab ~ Vaa) ~ (r]a ~ Vb)(Va ~ Vaab)} 

+ 2 Y^KU ~ 0(U - tbb) + (£« - 6) (6 - Ub)} 
P 

X [(Vab ~ r])(Vab ~ Vbb) + (Va ~ Vb){Vb ~ Vabb)], 



(61) 
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• the tadpole graphs: 
(III C ) = (/c a - k) 

X { (Cafe + Ccd ~ 2C) [ - (£a - £fe)(?7aa + ?7ad ~ Vbb ~ Vcb) - (Caa + £ad ~ £fefe ~ 6*)0?a ~ %)] 

P 

+ (Cad + Ccfe - 2C) [ - (£afe - O^aafe + Vabb ~ T} c - T)d) ~ (£aafe + £afeb ~ £c ~ £d) (%6 ~ 7?)] 

P 

+ (Ca - Cfe + Cc - Cd) 2J [(^ ~~ 0(^aa + Vad ~ Vbb ~ Vcb) + (£aa + £ad ~ £feb ~ £cb)(Vab ~ V) 
P 

+ (£a - ib){riaab + Vabb ~ % ~ Vd) + (£aafe + £abb ~ £c ~ £d)(Va ~ 

(62) 

• graphs with six external legs in the FT definition: 

(IV C ) = 4(K a + K ab -2K) [(C ab+ C cd -2C)(Cad+Ccfe-2C)-(Ca-Cfe+Cc-Cd) 2 ] , (63) 

• the BL difference: 



(V f ) = -A{2K a + K ab -?>K) (U + C«i-2C)(C««i + Cc6-2C)-(C«-C6 + Cc-Crf) 



(6- 



Note that (IV^) and (V^) do not involve sums over lattice points, and 
that they are proportional to the tree level expression (|58|). Their sum reads 

(iv c ) + (v ( ) = -4( Ka - K )[(c a fe+C cd -2C)(Ca d +C C fe-2C)-(Ca-Cfe+C c -0) 2 

(65) 

The simplicity of these expressions might be a hint of the existence of a 
more direct derivation. 



2.4 Numerical evaluation 

The procedure we employ to evaluate the correlator numerically follows the 
conventional philosophy adopted in formal perturbation theory: the contri- 
butions to each order are first evaluated in a finite volume and then the 
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Figure 4: Ranges for the computation: L determines the range of |x| for 
which topological charge correlators (? x are computed. L Q is a cutoff on the 
sum over Green's functions G x . L\ determines the range for which exact 
values of G x , as opposed to asymptotic approximations, are used. 

thermodynamic limit is taken order by order. It is hoped (though far from 
proven) that by this procedure one obtains an asymptotic expansion of the 
infinite volume correlator. For a more detailed discussion of the difficult 
issues involved see for instance jH]. 

The various perturbative contributions to the charge density correlator 
Go,x were computed, by means of C programs, for X in a finite 400 x 400 
sublattice of the whole 1R 2 . Lattice symmetry allows to reconstruct all corre- 
lators in the 2L x 2L square lattice (in our case, L = 200) from those in the 
shaded region. 

As explained in the next sections, the computation of Q x involves sums 
over Green's functions G x with |x| arbitrarily large, and a cutoff Lo on the 
summation range is imposed. Also, since we have an asymptotic expres- 
sion for G x , we only need to compute Green's functions exactly for |x| in a 
smaller region determined by another cutoff Lj; for larger |x| the asymptotic 
expression is used. 
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Figure 0] shows the ranges L, Li, Lq involved in the computation. 

2.4.1 Computation of Green's functions 

It is impractical to compute and store an array of all possible values of 
(subtracted) Green's functions G x according to the exact expression (pTT|). It 
is far more convenient to store only such values as differ noticeably from the 
approximate expression (|32j). We chose to compute exact values for G x with 
|x| less than a certain value Li (see figure EJ), and approximate values for the 
remaining ones as they were needed. 

The constant appearing in must be computed for this procedure to 
be useful. 

Working with L\ = 1000, exact Green's functions were computed 

for the triangular region < x < Li, < y < x. All those in the region 
|x| < L\ can be obtained from them by symmetry. 

Exact values of L\ at the edge x — L\ of the triangle were used to fit 
the Green's function space dependence to (jH2J), since this edge is the furthest 
removed from the centre and therefore provides the best agreement with the 
asymptotic form. 

A value c = 0.257343 is obtained for the constant (see figure EJ), which 
agrees with the analytical expression (J3l?|) taken from [S|. 

For Li = 1000, the numerical results for the correlators are stable with 
respect to 10% changes in L\. 

2.4.2 Computation of correlators 

The tree level correlator (jSHJ), as well as one- loop contributions (jo^j) and (jHU), 
do not involve sums over lattice points, and are readily computed for each 
X. 

One-loop contributions (joTIj) . (JoTj) and (|U2*|) involve sums over points arbi- 
trarily far away from this sublattice, and a cut-off must be introduced. We 
choose to restrict the sum to points within the circle of radius Lq about the 
origin. At Lo = 2000, results are stable with respect to 10% changes in Lq 
or in the shape of the summation domain (circle vs. square). This indicates 
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Figure 5: Fit of Green's functions to const - 7^:ln|x| for x = (Li, y), < 
y < L\ = 1000. The dashed line corresponding to the fitted dependence is 
indistinguishable from the computed points in this range. 
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the absence (or cancellation) of any infrared divergences connected with the 
summation over intermediate lattice points. 

Results are presented for the computations outlined in the previous sec- 
tions. Range parameters are L = 200, L\ = 1000, L Q = 2000. 

Figure El shows the (nonuniversal) near-origin tree level and one loop con- 
tributions to the topological charge density correlator Go,(x,o) as a function 
of x. Note that both contributions are positive at the origin, negative every- 
where else, and are suppressed for large x. The peak at the origin results, in 
the continuum limit, in the singular contact term discussed e.g. in JU|. The 
values of the peak at tree level and one loop level are 

Ql « 0.810569, Cfc « 0.899007. (66) 

The tree level peak can also be calculated analytically from ()3ip and 1)580 . 
yielding = 8/tt 2 « 0.810570. 

Since this is a perturbative analysis, the topological susceptibility (JHJ) 
should vanish order by order in (3~ x . Summing the computed correlators 
over the whole 2L x 2L region, the following results are obtained: 

w 1.3x 10~ 5 , J^S~3.7x 10" 5 , (67) 

X X 

which are compatible with zero and serve as a check of our computations. 
Note, in particular, that the problems deemed to arise with xt in the contin- 
uum limit of the 0(3) model do not show up in this perturbative treatment. 

Figure [7| shows the tree level contribution multiplied by function 
of x. This product approaches for large x the constant -0.2026, compatible 
with (jSnj), so the tree level result (including all numerical factors) is 

£o T (,,o) - -32^4 (68) 

for x 3> 1. 

The one-loop contributions (|50|l through (J64|) to the correlator were com- 
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Figure 6: Near-origin behavior of correlator Go,( x ,o)-> & t tree level and one 
loop, as a function of x. 
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Figure 7: Tree level Go,{x,o) multiplied by x 4 as a function of x. 



puted and their behavior for large x analyzed: 



(ic) 
(n c ) 
(III C ) 



x* 

blnx 



(IV C + V, 



a 

+ — > 

ar 



.T > 1, 



(69) 



with a = 0.101, b = 0.0653, d = 0.101. These Ansatze were obtained by 
inspection of the plots and trial and error, except for the last one, which 
stems directly from the proportionality of to (|58jl and the fact that 
k = 0, K a = — 1/4, which results in the asymptotics in ()69j) with d = 1/tt 2 . 

The leading terms in (1^) and (III^) cancel, but subleading terms 0(x~ 4 ) 
and 0(x~ A lnx) remain. The sum of all one- loop contributions (see figure |HJ) 
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Figure 8: One-loop level Go,(x$) multiplied by x A as a function of lnx, together 
with a linear fit for large x. 



behaves as 



Go A x 



c + Ci In x 



(*,o) 



x > 1, 



(70) 



with constants 

c = 2.206, ci = 0.0645. (71) 

Note that these are not identical to b and d in ()69j) . because the subleading 
terms in the asymptotic behavior of (If) and (III^) contribute. 

The same analysis was performed along the diagonal of points (x, x) with 
identical results for the asymptotic forms of the correlator: full rotational 
invariance is restored in this regime. 

2.4.3 Renormalization and continuum limit 



We now discuss the continuum limit, again in the spirit of formal perturbation 
theory, i.e. termwise in the perturbative expansion; no claim can be made 
that the resulting expansion is asymptotic to a nonperturbatively defined 
continuum limit. 
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Up to one loop we have found for the correlator on the unit lattice and 
large lattice distances |x| 3> 1 

Q - ~ ~w " W s o^W (co + Cl ln |x|) + (-f) ■ (72) 

To obtain a continuum limit we have to introduce the lattice spacing 
a, rescale x = y/a and make /3 dependent on the cutoff a. We also have to 
rescale the correlator Q by a factor a 4 according to its engineering dimension, 
to obtain the correlator in continuum normalization 

Q% = a 4 G*a ~ -— ^ — 1— --L- Ft4-^ U> + d In ^ ] i 



y xa /3(a) 2 327r 4 |y| 4 /3(a) 3 647r 2 |y| 4 V a / V/3( a ) 4 

(73) 

valid asymptotically for |y| ^> a. 

A continuum limit should exist if we let f3 depend logaritmically on a 
according to the one-loop Callan-Symanzik (3 function for the 0(3) model 

m = ^Xl^) + OV?) = 1 (l + ^ ln( H + O^)) . (74) 

Inserting this in ()73|) and reexpanding to order /J^ 3 we see that the terms 
proportional to lna cancel if c% = 2 /it 3 , which is consistent with the value 
0.0645 produced by our numerical computation; the continuum limit of the 
correlator is then 



\y\fie 2 



7T 3 £0 



+ O(/3 - 4 ) . (75) 

By choosing a renormalization scale \x = ^ we obtain 

lime; S 5y = --i ^ (l + -L >n Wy |)) + O (A -4 ) . (76) 

If we had chosen the FT instead of the BL definition, the only difference 
wuld have been a change in cq, corresponding to different renormalization 
scale, so to the order considered the two definitions are related by a finite 
renor malizat ion . 
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Equation (J76|) is now valid for all y, since in the limit a — > only the 
asymptotic behavior of (? x survives. 

We can compare this result with the Renormalization Group improved 
tree level result (21 El 



vy ~ 32n^(yy |y|« ' y J 

if we reexpand (J77j) to order /3q 3 , using the one-loop RG flow 

/3(y) = /?(yo)-^-lnM (78) 
^ |yol 

The expressions (f7U|) and (J77)l agree to order /3q 3 = /3(yo)~ 3 if the renormal- 
ization scales are chosen appropriately. 

So we have established that our one-loop calculation supports the soften- 
ing of the short distance behavior predicted by the RG improved tree level 
result. 



3 Conclusions 

We computed the two-point functions of the topological charge density per- 
turbatively to one loop and found consistency with the RG improved tree 
level perturbative result, indicating a softening of the short distance singu- 
larity compared to the naive tree level result. This means that at the level of 
formal one-loop perturbation theory the requirements of the two positivities 
analyzed in jT] are indeed satisfied. 

But we should point out again that a mathematical justification of the 
procedures of formal perturbation theory (interchange of the weak coupling 
limit with the thermodynamic and continuum limits) does not exist, there- 
fore it would be very interesting to see if this softening is really present at 
the nonperturbative level. Numerical checks using Monte Carlo simulations 
are presumably very difficult, however, since the check requires identifying 
logarithmic corrections to a rather large power. 
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